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Estimation  of  the  Range  between  the  Intercepter  and  the  Target 
during  Infrared  Terminal  Guidance  of  Space  Interception 


Shi  Xiaoping,  Wang  Zicai,  Ke  Qihong 

(Simulation  Center,  Harbin  Institute  of  Technology,  Harbin, 
150001) 


Abstract:  A  proper  mathematical  model  is  constructed 

according  to  the  kinematic  equations  of  space 
interception.  After  the  component  of  the  relative 
velocity  normal  to  the  line  of  sight  is  observed  using 
a  nonlinear  observer  design  method  with  two  canonical 
forms,  the  problem  of  estimating  the  range  between 
interceptor  and  the  target  during  infrared  terminal 
guidance  of  space  interception  is  solved.  Numerical 
simulation  verifies  the  above  results. 

Key  Words:  Space-interception,  terminal  guidance, 
optimal  control 


During  infrared  terminal  guidance  for  space  interception, 
the  relative  range  between  interceptor  and  target  cannot  be 
directly  measured.  Nonetheless,  it  has  to  be  estimated  since 
some  guidance  laws  virtually  rely  on  this  physical  quantity  [1]. 
As  far  as  nonlinear-observer  design  methods  are  concerned,  most 
are  subject  to  extremely  rigorous  application  conditions .  Here 
in  this  paper,  a  design  method  of  this  kind,  based  on  two 
canonical  forms,  is  introduced,  which  has  proved  to  be  applicable 
to  the  majority  of  nonlinear  systems.  Although  the  observer  that 
our  method  produces  may  possibly  contain  a  control  variable 
derivative,  a  closed-loop  observer  without  any  control  variable 
and  its  derivative  can  be  created  as  long  as  the  feedback  control 
law  is  a  continuous  differentiable  function  of  the  state.  The 
computing  load  of  the  foregoing  design  method  depends  on  the 
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order  and  structure  of  the  model.  An  appropriate  model  is 
presented  in  this  paper,  which  serves  in  designing  an  observer  to 
observe  the  relative  velocity  component  in  the  direction  normal 
to  the  line  of  sight,  thus  indirectly  solving  the  problem  of 
estimating  the  relative  range  between  interceptor  and  target. 


1 .  Constructing  a  Mathematical  Model 


The  relative  kinematic  equation  of  interceptor  and  target 
is  [1] 


r  =  V 

2 

V  =  rto 

2va)  “i 


(1) 


where  r  is  relative  range;  v  is  relative  velocity;  o  is  visual 
angular  velocity;  Uj  is  mobile  acceleration  of  the  interceptor. 
An  infrared  seeker  head  can  only  provide  w,  whose  value  is 
generally  at  the  10"*rad/s  order  of  magnitude.  In  addition,  v 
varies  slightly,  being  basically  a  constant  and  known.  From  Eg. 
(1),  the  following  can  be  derived; 


^=ra)’--(2va>  +  M,) 
dt  r  ' 


(2) 


and  since  state  and  control  variables  Xi=re,  X2=vtt, 

u=v(m>,  u=ve+Ui,  are  introduced.  As  a  result,  a  bivalent 
approximate  model  is  obtained  as  follows: 


=  -u=f^(x^,  x^,  u) 

•^2  ""  ~7"(^2  ■'■“)==/2(-*l*  ^2’  “) 

I 

y=x^  =h{x^,  Xy  u) 


(3) 
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where  the  measurement  equation  selected  is  rational.  According 
to  Eq.  (3),  the  observer  is  now  designed  to  observe  the  relative 
velocity  component  in  the  direction  normal  to  the  line  of  sight 
Xi,  and  then  to  estimate  r.  The  visibility  matrix  [2]  of  this 
system  is 

e(r,,  Xj,  «)  = 

Thereby,  the  systematic  equation  (3)  is  visible  as  long  as 


Xj+X^U 


2Xj  +M 


(4) 


2.  Design  of  Nonlinear  State  Observer 

For  the  convenience  of  design,  first  suppose  the  canonical 
form  I  of  the  systematic  equation  (3) 

X,‘  =Xj  =/’(x*  ,  Xj,  «,«.«) 

X*  =  —  a*(x*,  x‘,  u,  u,  u)=f^{x^  ,  Xj  .  II.  u,  u) 
y  =  h'  {xl ,  x’)  =  x’ 

The  conversion  from  Eq.  (3)  to  canonical  form  I  is^^^ 

9,(x,,  Xj,  II,  u)  =  h{x^y  x^,  ii)  =  x’ 

9j(x,,  Xj,  II,  u)  =  Nq^(x^,  ofj,  tt,  ii)  =  x* 

where  the  operator  N 

m 

lou  duJ 

where  f=[fi,  for  the  given  problem  in  this  paper,  the 

conversion  equation  ( 6  )  is 


+  riL 

L5x,  dxj^ 


(7) 
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xt  =  v,(x,,  Xj,  «.  «)|  (11) 

x^=Vj(x,,  Xj,  «)  i 

where  [x^,  X2]^6G,  and  Vj,  V2  are  continuous  differentiable 
functions  of  Xi,  X2,u  and  u  ;  the  total  derivative  with  regards  to 
t  in  the  foregoing  equation  should  be  solved  and  Eg.  (10)  should 
be  taken  into  consideration.  After  all  these  procedures,  the 
following  can  be  finally  arrived  at: 

0  =  Mc^(y,  u)  +  Ma^{y,  «,  u)  +  a^(y,  u,  li.  u) 

v,(x,,  X,,  M,  u)  =  Mc(y,  u)  +  a^(y,  u,  u) 

VjCx,.  Xj,  «)  =  Xj=c(>',  u) 

where  operator  M  is 

-Hi;  |][" 

where  C  is  a  scalar  continuous  differentiable  function. 

Combining  Eqs .  (12)  with  (5)  and  (9),  the  following  is  derived: 
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0  =  M*^c(y,  u)  +  M  a  (y,  u,  u)  +  a^{yf  u,  u,  ii) 


where  operator  M*  is 


,T 


where  f  “I^i  /j  1  ;  by  substituting  Eq.  (5)  in  the  foregoing 

equation  and  then  in  Eq.  (14),  the  following  is  acquired; 


,  A-;,  u,  «)  + 

2  2  Oy  I 


c(y>  w; 


0=  >  ^'2 


dydu  ^ 


(IS) 

da,{y,  u,  it)  da^(y,  u,u)  . 

— ! — 5 - u  +  — - u+aA}’,  u,  u,  u) 

Ou  Cu 

By  substituting  Eq.  (9)  in  the  foregoing  equation,  the  following 
partial  differential  coupled  equations  can  be  obtained  as 

follows;  ^ 

2  dc(y,  u)  _c~c(v’,  u) 

y  +  u  dy 

_ u  dc(y,  u)  ._2c'c(y.  u)  .  “) 

y  +  u  8y  dydu  ^  du 

.  ■\_dc(y,  u)  ,  ^«,Cv.  w.  u)  ■  . 

-a,0-,  «,  «,  u) - - y. - «  + 

da^(y,  u,u)  d\(y, 

Su  du-  -I 

From  the  foregoing  coupled  equations ,  a  set  of  simplest  solutions 
can  be  derived; 


c(y,  u) 


y  +  u 


a.{y,  u,  m)  = 


(f  +  «) 


fl„(v,  u,  li.  «)—  r 

Cv  +  w)  (y+w) 


2  ,..2 
— j(m) 
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Substituting  the  foregoing  equation  in  Eq.  (12),  the  conversion 
equation  (11)  is  acquired 


_ u 

2 

Xj  +  u 


(18) 


By  substituting  Eq.  (17)  in  Eq.  (10),  the  specific  canonical  form 
II  is  obtained  as  follows; 


.  A 

X,  =  — 


+  — ^(«)- 


'  (V  +  m)^  0  +  «) 


•  A  A 

Xj  =x,  + 


3« 


(y  +  uY 


-  y  +  u 


(19) 


The  observer  of  the  systematic  equation  (19)  is 


*A 

X  ,  =  — 


2  ,2  ,  -.A^ 

+ - r(M)  TT +^2) 


0  +  h)^  (y  +  uY 


y  +  u 


lA  ..A 
X,  =X,  + 


3ii 


(y  +  u)' 


(20) 


By  subtracting  Eq.  (20)  from  Eq.  (19),  the  state  error  equation 
is  obtained; 

;,A  ..A 

=  -Pl^2 

•  A 

X,  =x,  -P,x^ 

where  "  ’  °  The  parameters  selected  are  Pi=AiA2/ 

P2~“(Ai+A2),  where  Aj  and  A2  are  poles  of  the  systematic-error 


(21) 
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equation  (21).  Under  the  action  of  the  conversion  equation  (18), 
the  observer  of  systematic  equation  (3)  can  be  acquired  from 
Eq.  (20) 


ix^+u)x 


X.  =■ 


2i, 


.\J(^ _ — 

Lcv  +  u)^  {y  +  uf 


{  1  . 

^  'll 

Vxj+M 

y +  m/ 

(u'  +  MX,  —  MX|)x, 

lE 


3m 


2Jc, 


(y  +  M)  'V,(x2+k)  (x, +m) 


+ 


+  - 


(ux^  —  2mXj  —  2mm)x, 


2Xj(Xj  +«) 


(x,  +m) 


2  r- 


^2=- 


3m 


2^2 


Lo  +  m)'  x,C^j  +  m)- 


- ^)]-h 

^VXj+M  y  +  M/J  2 


(22) 


3.  Estimation  of  Relative  Range  between  Interceptor  and  Target 


To  eliminate  m>'  «'  and  u  in  the  foregoing  equation,  the 
yisual  equation  (3)  is  a  monovalent  system.  In  accordance  with 
the  principle  of  designing  a  guidance  law,  r  and  e  should  be 
made  to  approach  zero  so  as  to  give  the  index 


1  r  V 

2J  [9,^!(o+j,«^(0]<i^ 


(23) 


then  the  optimal  control  can  be  solved  as 

ui.t)  =  yPi.f)xS^) 


(24) 


where  p(t)  can  satisfy  Riccati  differential  equation 


(25) 
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The  solution  to  that 


the  final  value  condition  is  p(tf)=ri. 
equation  is 


p{t)  = 


(26) 


Based  on  Eqs.  (3),  (6)  and  (25),  the  following  can  be 
derived : 

«(0=  r  (/) 

s. 


(27) 

(28) 


From  Eqs.  (24),  (27)  and  (28),  select 


u  =  ^p{t)x^ 


u  =^■p{t)x^ 
5. 


(29) 


By  substituting  the  foregoing  equation  in  Eq.  (22),  the  closed- 

loop  observer  can  be  constructed.  The  visual  relative  velocity  v 
is  a  constant,  then  from 


Xj 


=  fo) 


(30) 


the  relative  range  can  be  estimated  as 


O)  VO)  ^  . 


Fig.  1  shows  the  structure  of  state  estimation  and  the  generation 
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of  relative  range  estimation. 


4.  Numerical  Simulation  Check 

During  a  numerical  simulation,  select  ri=0.4,  qi=0.3,  Si=0.3, 
Ai=-20,  X2=-30.  Selecting  to=0,  tf  is  the  blind  region  entry- 
moment,  r(0)=150km,  v=-7.5km/s,  ©( 0 )=5xl0‘*rad/s ,  then 
Xi(0)=75m/s,  X2( 0 )=-3 . 75m/s^;  select  x,(0)=80m/s,  jCj(0)=  —  4m/s\ 
Figures  2-5  show  the  desirable  numerical  simulation  results. 


1 


Fig.  1.  Structure  of  State  Estimation 
Key:  1.  Model  D;  2.  State  observer 


Fig.  2  State  and  Its  Estimated  Value 
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T/s  T/s 


Fig.  3  State  X2  and  Its  Estimated  Value 


Fig.  4  Relative  Range  r  and  Its  Estimated  Value  f 


Fig.  5  Visual  Angular  Velocity 


5 .  Conclusions 

By  using  our  nonlinear  state  observer  design  method  in  two 
canonical  forms,  estimation  of  the  relative  range  between 
interceptor  and  target  during  infrared  terminal  interception  was 
indirectly  made  through  an  appropriate  model  and  confirmed  by 
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numerical  simulation  results.  As  for  the  stability  of  the 
closed-loop  nonlinear  feedback  system  with  an  observer,  see 
reference  [3]. 
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